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Abstract 

We calculate the electric charge radii of the SU (3) pseudoscalar mesons and the SU (3) octet 
baryons in quenched and partially quenched chiral perturbation theory. We work in the isospin 
limit, up to next-to-leading order in the chiral expansion, and to leading order in the heavy baryon 
expansion. The results are necessary for the extrapolation of future lattice calculations of meson 
and baryon charge radii. We also derive expressions for the nucleon and pion charge radii in SU(2) 
flavor away from the isospin limit. 
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I. INTRODUCTION 



The study of hadronic electromagnetic form factors at low momentum transfer provides 
important insight into the non-perturbative structure of QCD. A model- independent tool to 
study QCD at low energies is chiral perturbation theory (%PT), which is an effective field 
theory with low-energy degrees of freedom, e.g., the meson octet in SU(3) flavor or the pion 
triplet in the case of SU(2) flavor. assumes that these mesons are the pseudo-Goldstone 
bosons that appear from spontaneous breaking of chiral symmetry from SU(3)l <8> SU(3)r 
down to SU(3)v- Observables receive contributions from both long-range and short range 
physics; in ;^PT the long-range contribution arises from the (non-analytic) structure of pion 
loop contributions, while the short-range contribution is encoded in a number of low-energy 
constants that appear in the chiral Lagrangian and are unconstrained in %PT. These low- 
energy constants must be determined from experiment or lattice simulations. 

Notable progress toward measuring the proton and neutron form factors has been made 
in recent years and high precision data are available (see 0,0] for references). Experimental 
study of the remaining octet baryons, however, is much harder. The charge radius of the 
X~ has only recently been measured Although more experimental data for the other 
baryon electromagnetic observables can be expected in the future, progress will be slow 
as the experimental difficulties are significant. Theory, however, may have a chance to 
catch up. While quenched lattice calculations have already appeared 0, IE IS 0-IEE1- with 
the advance of lattice gauge theory, we expect partially quenched calculations for many of 
these observables in the near future. One problem that currently and foreseeably plagues 
these lattice calculations is that they cannot be performed with the physical masses of 
the light quarks. Therefore, to make physical predictions, it is necessary to extrapolate 
from the heavier quark masses used on the lattice (currently on the order of the strange 
quark mass) down to the physical light quark masses. For lattice calculations that use the 
quenched approximation of QCD (QQCD), where the fermion determinant that arises from 
the path integral is set to one, quenched chiral perturbation theory (Q%PT) [13, EH E2, EH 
3 3 3] has been developed to aid in the extrapolation. The problem with the quenched 



approximation is that the Goldstone boson singlet is no longer affected by the U(1)a anomaly 
as in QCD. In other words, the QQCD equivalent of the rf that is heavy in QCD remains light 
and must be included in the Qx?T Lagrangian. This requires the addition of new operators 
and hence new low-energy constants. In general, the low-energy constants appearing in the 
QxPT Lagrangian are unrelated to those in %PT and extrapolated quenched lattice data is 
unrelated to QCD. In fact, several examples show that the behavior of meson loops near the 
chiral limit is frequently misrepresented in QxPT HI, HI 13] • We find this is additionally 
true for the meson and baryon charge radii. 

These problems can be remedied by using partially quenched lattice QCD (PQQCD). 
Unlike QQCD, where the masses of quarks not connected to external sources are set to 
infinity, these "sea quark" masses are kept finite in PQQCD. The masses of the sea quarks 
can be varied independently of the valence quark masses; usually they are chosen to be 
heavier. Sea quarks are thereby kept as dynamical degrees of freedom and the fermion 
determinant is no longer equal to one. By efficaciously giving the sea quarks larger masses it 
is much less costly to calculate observables than in ordi nary QCD. The low-energy effective 
theory of PQQCD is PQxPT 0, HI 0. Since PQQCD retains a U(1) A 
anomaly, the equivalent to the singlet field in QCD is heavy (on the order of the chiral 
symmetry breaking scale A x ) and can be integrated out (2ll. U^. Therefore, the low-energy 
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constants appearing in PQ^PT are the same as those appearing in xPT. By fitting PQ^PT 
to partially quenched lattice data one can determine these constants and actually make 
physical predictions for QCD. PQyPT has been used recently to study heavy meson 29] 
and octet baryon observables |3(1|31> 32l . 

While there_are a number of lattice calculations for observables such as the pion form 

36l | that use the quenched 



factor |4j, |33j, |34j or the octet baryon magnetic moments [3J 
approximation, there are currently no partially quenched simulations. However, given the 
recent progress that lattice gauge theory has made in the one-hadron sector and the prospect 



of simulations in the two-hadron sector 37, |H, |3^, H(l 41 1, we expect to see partially 



quenched calculations of the electromagnetic form factors in the near future. 

The paper is organized as follows. First, in Section |HJ we review PQxPT including 
the treatment of the baryon octet and decuplet in the heavy baryon approximation. In 
Section ITTT1 we calculate the charge radii of the meson and baryon octets in both QxPT and 
PQxPT up to next-to-leading (NLO) order in the chiral expansion. We use the heavy baryon 
formalism of Jenkins and Manohar 0, HH , treat the decuplet baryons as dynamical degrees 
of freedom, and keep contributions to lowest order in the heavy baryon mass, Mb- These 
calculations are done in the isospin limit of SU(3) flavor. For completeness we also provide 
the PQxPT result for the charge radii for the SU(2) chiral Lagrangian with non-degenerate 
quarks in the Appendix. In Section |IV| we conclude. 



II. PQxPT 

In PQQCD the quark part of the Lagrangian is written as 0, [i^, HH, BE SHE 

£ = ^(ilp- m q ) ab q b + ^ Qa{ i P- m q)ah%+ &ea, a (^-ra S eaW Se a,6 

a,b=u,d,s a,b=u,d,s a,b=j,l,r 

= ^2 Qj( iJ P- m o)jkQk- (i) 

j ,k=u,d,s ,u,d,s ,j ,l,r 

Here, in addition to the fermionic light valence quarks u, d, and s and their bosonic coun- 
terparts u, d, and s, three light fermionic sea quarks j, I, and r have been added. These 
nine quarks are in the fundamental representation of the graded group 
and have been accommodated in the nine-component vector 

Q = (u,d,s,j,l,r,u,d,s) (2) 

that obeys the graded equal-time commutation relation 

Q?(x)Qf (y) - (-i)**qf (y)Of (x) = rt^ 3 (x - y ), (3) 

where a and (3 are spin and % and j are flavor indices. The graded equal-time commutation 
relations for two Q's and two Q^s can be written analogously. The grading factor 

f 1 for k = 1,2,3,4,5,6 , . 

Vk ~ \ for k = 7, 8, 9 W 

takes into account the different statistics for fermionic and bosonic quarks. The quark mass 
matrix is given by 

m Q = diag(m„, m d , m s , rrij, mi, m r , m u , m d , m s ) (5) 
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so that diagrams with closed ghost quark loops cancel those with valence quarks. Effects 
of virtual quark loops are, however, present due to the contribution of the finite-mass sea 
quarks. In the limit rrij = m u , mi = rrid, and m r = m s QCD is recovered. 

It has been recently realized {47] that the light quark electric charge matrix Q is not 
uniquely defined in PQQCD. The only constraint one imposes is for the charge matrix Q 
to have vanishing supertrace. Thus as in QCD, no new operators involving the singlet 
component are subsequently introduced. Following (3(| we use 

Q = diag q-j, q h q r , q h q h q^j (6) 

so that QCD is recovered in the limit mj — > m u , mi — > m^ and m r — > m s independently of 
the g's. 



A. Mesons 

For massless quarks, the Lagrangian in Eq. exhibits a graded symmetry SU(6\3)l <S> 
SU(6\3)r Cg> U{l)y that is assumed to be spontaneously broken down to SU(6\3)y <%> U{l)y. 
The low-energy effective theory of PQQCD that emerges by expanding about the physical 
vacuum state is PQxPT. The dynamics of the emerging 80 pseudo-Goldstone mesons can 
be described at lowest order in the chiral expansion by the 0(E 2 ) Lagrangian 1 

C = J —stx (D^D^) + A str (m Q £ + m^) + ad^ d^ - fifol (7) 

where 

S = exp(^) =?2 , (8) 



X M 

f = 132 MeV, and we have defined the gauge-covariant derivative D^T, = (9 M S + ieA^[Q, £]. 
The str() denotes a supertrace over flavor indices. The M, M, and x are matrices of 
pseudo-Goldstone bosons with quantum numbers of qq pairs, pseudo-Goldstone bosons with 
quantum numbers of qq pairs, and pseudo-Goldstone fermions with quantum numbers of qq 
pairs, respectively. $ is defined in the quark basis and normalized such that $12 = vr + (see, 
for example, Upon expanding the Lagrangian in (JZj) one finds that to lowest order the 

mesons with quark content QQ' are canonically normalized when their masses are given by 

4A 

m QQ' = -p( m Q + m Q')- ( 10 ) 

The flavor singlet field given by $ = str($)/ v / 6 is, in contrast to the Qx?T case, rendered 
heavy by the U(1)a anomaly and can therefore be integrated out in %PT. Analogously its 
mass /io can be taken to be on the order of the chiral symmetry breaking scale, fiQ — > A x . In 
this limit the flavor singlet propagator becomes independent of the coupling a and deviates 
from a simple pole form |2l|, |22j . 



Here, E ~ p, where p is an external momentum. 
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B. Baryons 



Just as there are mesons in PQQCD with quark content QiQj that contain valence, 
sea, and ghost quarks, there are baryons with quark compositions QiQjQk that contain all 
three types of quarks. Restrictions on the baryon fields Byk come from the fact that these 
fields must reproduce the familiar octet and decuplet baryons when i, j, k = 1-3 [3, l3fll |. 
To this end, one decomposes the irreducible representations of SU(6\3)y into irreducible 
representations of S77(3) va i <8> S77(3) sea <8> S7^(3) g host <8> U(l). The method to construct the 
octet baryons is to use the interpolating field 

Bl jk - (QT a Qf b Ql' c - Q? a Q] ,c Qk h ) ^(c^U (n) 

which when restricted to i, j, k — 1-3 has non-zero overlap with the octet baryons. Using 
the commutation relations in Eq. Q one sees that B^k satisfies the symmetries 

Bijk ( ) Bikji 

= B ijk + (-) l+r > iT >*B jik + (-y+^+mv k +v k vi Bkj ._ (12 ) 

The spin- 1/2 baryon octet B^f. = B^k, where the indices i, j, and k are restricted to 1-3, 
is contained as a (8, 1, 1) of SU (3) va i ® SU (3) sea ® 5?7(3)ghost in the 240 representation. The 
octet baryons, written in the familiar two-index notation 

/76 A + 75 S ° S+ 
B=l S- ^A-^ S o n |, (13) 

are embedded in as [H| 

-Bijfc = ( e ijiBki + e^iBji) . (14) 

Besides the conventional octet baryons that contain valence quarks, qqq, there are also 
baryon fields with sea and ghost quarks contained in the 240, e.g., qq se a(j- Since we are only 
interested in calculating one-loop diagrams that have octet baryons in the external states, 
we will need only the B^k of two valence and one sea quark or two valence and one ghost 
quark. We use the construction in [30]. 

Similarly, the familiar spin-3/2 decuplet baryons are embedded in the 165. Here, one 
uses the interpolating field 

%7 - {QT a Qf b QT + Qi' b Q] ,c Qk a + QT c QT a Qk b ) e abc (cy% 7 (is) 

that describes the 165 dimensional representation of SU(6\3)y and has non-zero overlap 
with the decuplet baryons when the indices are restricted to i, j, k = 1-3. Due to the 
commutation relations in Eq. (J3*|). %jk satisfies the symmetries 

% jk = {-y~ im Tjik = n 1+wk ^kj. (i6) 

The decuplet baryons are then readily embedded in T by construction: = 7^, 

where the indices i, j, k are restricted to 1-3. They transform as a (10, 1, 1) under 
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S£/(3)vai <g> SU(3) sea <g> 5C/(3) ghost . Because of Eq. JEJ, Tjjfc is a totally symmetric ten- 
sor. Our normalization convention is such that Tm = A ++ . For the spin-3/2 baryons 
consisting of two valence and one ghost quark or two valence and one sea quark, we use the 
states constructed in [3(J. 

At leading order in the heavy baryon expansion, the free Lagrangian for the and %j k 
is given by 2| 

C = i(Bv- VB) + 2a M (BBM+) + 2(3 M (BM+B) + 2a M (BB) str (M+) 

-i (T% • V%) + A (T%) + 2 7M (T^M + %) - 2a M (T%) str (M+) , (17) 

where A4 + = | {^ m Q^ + C m Q^)- The brackets in (|17j) are shorthands for field bilinear 
invariants originally employed in [15]. To lowest order in the chiral expansion, Eq. (|17|) 
gives the propagators 

vk" vk-A { } 

for the spin-1/2 and spin-3/2 baryons, respectively. Here, v is the velocity and k the residual 
momentum of the heavy baryon which are related to the momentum p by p = Mbv + k. 
Mb denotes the (degenerate) mass of the octet baryons and A the decuplet-baryon mass 
splitting. The polarization tensor 

pv = ( v "v v - g^) - -S»S U (19) 

3 

reflects the fact that the Rarita-Schwinger field {T^ijk contains both spin-1/2 and spin-3/2 
pieces; only the latter remain as propagating degrees of freedom (see |42| . for example). 

The Lagrangian describing the relevant interactions of the B^k and %j k with the pseudo- 
Goldstone mesons is 



C = 2a (BS^BA^) + 2(3 (BS^A^B) + J|c [(TA V B) + h.c] . 



(20) 



The axial- vector and vector meson fields A 11 and are defined by analogy to those in QCD: 

and V" = i(^t + £t3^). (21) 

The latter appears in Eq. (j2*Uj) in the covariant derivatives of B^ and %jk that both have 
the form 

{V»B) ijk = ^ fc + (^i%, + (-)^^^ (22) 

The constants a and (3 are easily calculated in terms of the constants D and F that are 
used for the SU(3) va \ analogs of these terms in QCD. Restricting the indices of Bijk to 
i, j, k = 1, 2, 3 one easily identifies 

a = -D + 2F and (3 = --D + F. (23) 
3 3 
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FIG. 1: Loop diagrams contributing to the octet meson charge radii in PQxPT. Octet mesons are 
denoted by a dashed line, singlets (hairpins) by a crossed dashed line, and the photon by a wiggly 
line. Only the third diagram has q 2 dependence and therefore contributes to the charge radius. 



III. CHARGE RADII 



In this Section we calculate the charge radii in PQ^PT and Q%PT. The basic conventions 
and notations for the mesons and baryons in PQxPT have been laid forth in the last section; 
QxPT has been extensively reviewed in the literature 3 2 3|- 2 



A. Octet Meson Charge Radii 

The electromagnetic form factor Gx of an octet meson (fix is required by Lorentz invari- 
ance and gauge invariance to have the form 

(Mp')\JL\Mp)) = G x (q 2 )(p + P r (24) 

where g M = (p' — pY and p (p') is the momentum of the incoming (outgoing) meson. Conser- 
vation of electric charge protects it from renormalization, hence at zero momentum transfer 
eGx(0) = Qx, where Qx is the charge of (fix- The charge radius r x is related to the slope 
of Gx(q 2 ) at q 2 = 0, namely 

<r 2 x >=Q^- 2 G x (q 2 )\ q i =0 . (25) 

There are three terms in the 0(E 4 ) Lagrangian 

8A 8A 
C = a 4 ^str(^ t S^S)str(m Q S + m^S t ) +a 5 -^str( J D At S^S(m Q S + mJ ? S t )) 

+ia^ti{L^D^D v T) + R„ V D^D V T) + ... (26) 

that contribute to meson form factors at tree level. Here L^, are the field-strength 
tensors of the external sources, which for an electromagnetic source are given by 

V = = eQ{d^A u - d u A„) + ie 2 Q 2 [A ll ,A u ]. (27) 



2 In short, there are only valence and ghost quarks in the theory and they have the same mass and charge 
pairwise so that diagrams with disconnected quark loops are zero. The singlet is light and has to be 
retained in the theory. It is treated perturbatively and resulting observables generally depend upon the 
parameters a and /zq. 
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FIG. 2: Wavefunction renormalization diagrams in PQxPT. These diagrams, together with the 
third one in Fig. (^Q), ensure meson electric charge non-renormalization. 



Unlike Q%PT, where the low-energy constants are unique and have no known connection to 
%PT, in PQ^PT the parameters in ()26|) are the dimensionless Gasser-Leutwyler coefficients 
of xPT (4^] which can be seen by looking at mesons that contain sea quarks only. 

To calculate the charge radii to lowest order in the chiral expansion one has to include 
operators of i2 in (J7J) to one-loop order [see Figs. (JTJ) and (J2J)] and operators of (pjj) to tree 
level. Using dimensional regularization, where we have subtracted - + 1 — 7 + log47r, we 
find in PQxPT for the 7r + 



^- 2 [2F uj + F ur ]+a,-q\ 



(28) 



which interestingly does not depend on the charges of the sea and ghost quarks, qj, qi, q r 
For the K + we find 



1 + 



16tt 2 / 2 



- Qji )F UU - [-- q fl F % 



"j 



q r F ur + 



5 - Qji + q r ) F us - [- + qji } F js - - + q r ) F r 



+ a 9 j^q 2 (29) 



and for the K° we find 



G P K % 2 ) 



16tt 2 / 2 



Qji ) F uu + - + qji F uj + \ - + q r ) F ur + [ - + q r ) F. 



qji + q r ] F us - - + q^ F js - 



q r F r , 



(30) 



Here q^ = qj + qi and we have defined 



q m hn< o 



QQ' 



(31) 



where the function ^(a) is given by 



a 2 
6 ~ 3 



1 log 

a 



te - 1 



1 - ! + ie + 1 



5a 4 
f 18 ~ 3' 



(32) 



The first derivative of Fqqi at q 2 = 0, needed to calculate the charge radii, becomes 



r d t? 1 1 m QQ' , 1 



(33) 



Charge conjugation implies 



G*? = -G%, G£?=-G£?, and G$ = -G$, (34) 

which we have also verified at one-loop order. The form factors of the flavor diagonal mesons 
are zero by charge conjugation invariance. In the limit rrij — * fh, m r — > m s we recover the 
QCD result |48j, |49| as expected. 

It is interesting to note, that duplicating these calculations for Q%PT shows that there 
is no meson mass dependence at this order. Specifically we find 

(g 2 ) = (g 2 ) = G Q K+ (g 2 ) = -Gg_ (g 2 ) = 1 + ±a 9 V, (35) 

and the form factors of the neutral mesons are zero. Here we annotate the quenched constant 
a® with a "Q" since its numerical value is different from the one in Eq. (|26[). Eq. ()35|) reflects 
that flavor- singlet loops do not contribute to the g 2 -dependence at this order; thus the virtual 
quark loops are completely removed by their ghostly counterparts. This can readily be seen 
by considering the quenched limit of Eqs. (}28|) - (j30|) . The meson mass independence reveals 
once again the pathologic nature of the quenched approximation and seriously puts into 
question ;\PT extrapolations to the physical pion mass. 



B. Octet Baryon Charge Radii 

The electromagnetic form factors at or near zero momentum transfer that enable the 
extraction of the baryon ma gne tic moments and charge radii have been frequently inves- 
tigated in QCD ji^, IHO, 51. I52I E3 . I54I IHii \E$L 57, 58[. There are also recent quenched 
and part ially qu enched calculations of the octet baryon magnetic moments in QxPT and 
PQxPT [T9T 130. 159]]. Here, we extend these calculations to the octet baryon charge radii. 
We retain spin-3/2 baryons in intermediate states since formally A ~ m^. 

Using the heavy baryon formalism 0, |43| , the baryon matrix element of the electromag- 
netic current J^ 1 can be parametrized in terms of the Dirac and Pauli form factors F\ and 
F 2 , respectively, as 



(B{p') I J M | B{p)) = u{p') < v^F^q 



-qvF 2 (q 2 ) }u(p) 



with q = p' 



M B 

p. The Sachs electric and magnetic form factors defined as 

-,2 



(36) 



G E (q 2 ) 
G M {q 2 ) 



Fi(q 2 ) + 



q 



4M|" 



-') 



F 1 (q 2 )+F 2 (q 2 



(37) 
(38) 



are particularly useful. The baryon charge Q, electric charge radius < r\ >, and magnetic 
moment /i can be defined in terms of these form factors by 



< r 2 E >'- 



d 



Q = G E (0) 
Here the baryon charge Q is in units of e. 



6 df GEiq 



and fi = Gm(0) — Q. 



(39) 
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1. PQxPT 



Let us first consider the calculation of the octet baryon charge radii in PQxPT. Here, the 
leading tree-level correction to the magnetic moments come from the dimension-5 operators 3 

ip 

Z1V±b 

which can be matched on the QCD Lagrangian upon restricting the baryon field indices to 
1-3 

C= 2M~ B [^tr(B[S„S u ]{Q,B}) +f i F ti(B[S„S u }[Q,B])}F^ (41) 

where 

2 5 
A*a = + 2/iF and = + Hf (42) 

at tree level. The magnetic moments contribute the so-called Foldy term to charge radii via 
i^(0) in Eq. (}37j) . Likewise, further leading tree-level corrections to the charge radii come 
from the dimension-6 operators 

C = [c a (BBQ) + C/3 (BQB)] v,d v F^ (43) 

x 

and the parameters c+ and c_, defined by 

2 5 

c a = ~c + + 2c_ and cp = — c + + c_, (44) 

3 3 

are the same as those used in QCD. Here, we take the chiral symmetry breaking scale 
A x ~ Airf for the purpose of power counting. The NLO contributions arise from the one- 
loop diagrams shown in Figs. © and (@J). To calculate the charge radii we need the form 
factors Fi to first order in q 2 and F 2 (0) we find 

6 ,„ ,3 



<r E > = - — (Qc- + a D c + ) + 77772 (Q^ + a DfJ>D) 



A 2 ^ - " 2M 2 



16n 2 f 



x 



A x log ^ - 5 Px log + 10 (3' x G{m x , A, n) 



(45) 



Here, we have defined the function G(m, A,//) by 



m 



2 



A , A - VA 2 - m 2 + ie 



g(m, A, (j.) = log— - - = log- = . (46) 

^ V A 2 - m 2 A + V A 2 -m 2 + ie 

Note that in Eq. (J43j) the only loop contributions we keep are those non- analytic in m x . 

The parameters for the tree-level diagrams are listed in Table U] The computed values 
for the j3 x , f3' X) and A x coefficients that appear in Eq. (J43j) are listed for the octet baryons 
in Tables ITTl IIXI The corresponding values for the AS transition are given in Table IXl 
In each table we have listed the values corresponding to the loop meson that has mass m x . 
If a particular meson is not listed then the values for (3 X , j3' x , and A x are zero. 4 



3 Here we use = d^A v - dyA^, 

4 We have defined the coefficients (3 X and fix' to correspond to those defined in [30| where \x = Q \xf 
aj) /i£> + -j^rp J2x [PxnT-x + l3' x T{m x , A, /i)] and the function T{m Xl A, fx) is given in [3ol | . 
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FIG. 3: Loop diagrams contributing to the baryon magnetic moments and charge radii. A thin 
(thick) solid line denotes an octet (decuplet) baryon. The last two diagrams in the first row 
contribute to both magnetic moments and charge radii; the magnetic moment part of which has 
already been calculated in j^]. The first diagram in row 1 contributes q 2 dependence only to F\ 
and therefore is relevant for the charge radii. The remaining diagrams have no q 2 dependence. 
These along with the wave function renormalization diagrams in Fig. (JIJ maintain charge non- 
renormalization. 




X 



FIG. 4: Wave function renormalization diagrams needed to maintain baryon electric charge non- 
renormalization. 



TABLE I: Tree-level contributions in QCD, QQCD, and PQQCD 





Q 


CUD 


p, s+ 


l 


l 

3 


n, 3° 





2 
3 


S° 





1 
3 


£- 3- 


-1 


1 

3 


A 





1 
3 


S°A 





1 
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TABLE II: The coefficients j3 x , P' x , and A x in SU(3) flavor PQxPT for the proton. 



X 


0x 


0x 




7T 

K 

uj 
ur 


-\ (7D 2 + 6DF - 9F 2 ) - ± (577 2 - 6DF + 9F 2 ) g iZ 
-± (57> 2 - 677F + 9F 2 ) (1 + 3g r ) 
-| (D + 3F) 2 + | (5D 2 - 6DF + 9F 2 ) q jt 
-± (77 + 3F) 2 + ± (57> 2 - 677 F + 9F 2 ) q r 


(-1 + foO c 2 

(l8 + 6*0 C 2 

-±fl C 2 


-1 + 

1 + 3<? r 

2 - 3qji 
1 — 3qv 



TABLE III: The coefficients (3 X , /3' x , and A x in ST7(3) flavor PQxPT for the neutron. 



X 


(3 X 


0'x 


A x 


7T 

K 

uj 
ur 


i (1777 2 - 67>F + 9F 2 ) - | (57> 2 - 6/7F + 9F 2 ) g iZ 
-i (577 2 - 677 F + 9F 2 ) (1 + 3q r ) 
-§ (77 2 - 3/7F) + | (57> 2 - 677F + 9F 2 ) q jt 
-§ (£> 2 - 3/7F) + i (57> 2 - 677F + 9F 2 ) q r 


(h + M<? 
(fs + M c 2 2 
(§ - hji) c 2 

(l8 ~~ 6*0 ^ 2 


-1 + 3^z 
1 + 3q r 

-Ma 

-3q r 



2. Q X PT 



The calculation of the charge radii can be easily executed for Q%PT. The operators in 
Eqs. ()41|) and ()43j) contribute, however, their low-energy coefficients cannot be matched onto 
QCD. Therefore we annotate them with a "Q". Additional terms involving hairpins [lH Hit] 



do not contribute as their contribution to the charge radii is of the form (//q/A 
therefore of higher order in the chiral expansion. We find 



logm„ and 



<r%> 



(47) 



As with the meson case, the diagram where the photon couples to the closed meson loop does 
not contribute to baryon charge radii in the quenched case, cf., A x is zero. The remaining 
coefficients appearing in Eq. ()47|) are listed in Table IXII and stem from meson loops formed 
solely from valence quarks. 



TABLE IV: The coefficients fix, fix, and A x in 517(3) flavor PQxPT for the E+. 

fix fix A x 

I {D 2 + 3F 2 ) (1 - 3^) (-i + ± Qjl ) C 2 -| + 2q jt 

-i (lLD 2 + 6DF + 3F 2 ) - (D - F) 2 qjl - § (D 2 + 3F 2 ) q r + \q fl + ±q r ) C 2 \ + Qjl + 2q r 

-l{D-F) 2 {l + 3q r ) {^f + hr)C 2 \ + q r 

-| {D 2 + 3F 2 ) (4 - 3 Qjl ) - C 2 | - 2q jl 

-l (D 2 + 3F 2 )(2-3q r ) (^ - ±q r ) C 2 § - 2<? r 

l(D-F) 2 (2 + 3 qjl ) (_^_i^ 7 ) C 2 _2_ qji 

l(D-F) 2 (l + 3q r ) (-i- ±c/ r )C 2 



-.2 1 

3 
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TABLE V: The coefficients X , P' x , and A x in SU{3) flavor PQ^PT for the E°. 

0x Px A x 

§(£> 2 + 3F 2 )(l-3g j7 ) (_^ + ^)c 2 ~l + 2q 3l 

-\ (5D 2 + 6DF + 3F 2 ) — (D — F) 2 Qjl - § (D 2 + 3F 2 ) g r (-^ + §^ + ^g r ) C 2 \ + Qjl + 2q r 

-l(D-F) 2 (l + 3q r ) (^+i gr )C 2 \ + q r 

-| (l) 2 + 3F 2 ) (1 - 3 qjl ) (£ - ^-,) C 2 | - 2^-, 

-i(Z? 2 + 3F 2 )(l-6(7 r ) (_i__^ gr ) C 2 i_ 2(fr 

\{D-F) 2 {2 + 3 qjl ) (_^_i^)c 2 -1-^ 

i( J D-F) 2 (l + 3g r ) (_^_i gr )C 2 -\-q r 



TABLE VI: The coefficients (3 X , P' x , and A x in SU(3) flavor PQ^PT for the E". 

& ^ 

I + 3^ 2 ) (1 - 3^-,) (-£ + C 2 -| + 2^, 

i (D 2 - 6DF - 3F 2 ) - (D - F) 2 ^ - § (D 2 + 3F 2 ) g r (-£ + foj + ^g r ) C 2 § + q jt + 2q r 

-l(D-F) 2 (l + 3q r ) ^ + l Qr )C 2 \ + q r 

| (D 2 + 3F 2 ) (2 + 3fc-,) (-£ - ±q jt ) C 2 -f - 2g, 7 

|(D 2 + 3F 2 )(l + 3g r ) {-i.-y r )C 2 -l-2q r 

i(D-F) 2 (2 + 3 %! ) (_^_1^ 7 )C 2 -§-g j7 

i (jD _F) 2 (l + 3g r ) (_^_i gr )C 2 -I-g r 



IV. CONCLUSIONS 

We have calculated the charge radii for the octet mesons and baryons in the isospin limit 
of PQxPT and also derive the result for the nucleon doublet and pion triplet away from the 
isospin limit for the SU(2) chiral Lagrangian. For the octet mesons and baryons we have 
also calculated the QxPT results. 

Knowledge of the low-energy behavior of QQCD and PQQCD is crucial to properly 
extrapolate lattice calculations from the quark masses used on the lattice to those in nature. 
For the quenched approximation, where the quark determinant is not evaluated, one uses 



TABLE VII: The coefficients p x , P' x , and A x in SU(3) flavor PQ^PT for the S°. 



Px P x A x 



l(D-F) 2 (l-S qjl ) (_^+i^) C 2 _i +qji 

i (11D 2 + 6DF + 3F 2 ) -\{D 2 + 3F 2 ) Qjl - (D - F) 2 q r + ± qjl + \q r ) C 2 -\+q r + 2 Qjl 



r 



l(D 2 + 3F 2 ){l+3q r ) (A_ + A_ q J C 2 \+2q, 

-\ {D -F) 2 (A-3 q]l ) (4f-hn)C 2 i-Qji 



72 2 -a 

3 " r 



-l(D-F) 2 (2-3q r ) (wr-hr)C 2 
| (D 2 + 3F 2 ) (2 + 3 Qjl ) (-£ - C 2 -I - 2q 3l 

| (D 2 + 3F 2 ) (1 + 3q r ) (-^-^ qr )C 2 -\-2q r 
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TABLE VIII: The coefficients p x , P' x , and A x in SU{3) flavor PQ^PT for the E~. 



X 




Px 


Px 


A x 


TT 

K 

uj 
ur 

sj 
sr 


~\{D 2 


\{D-F) 2 {l + 3 qjl ) 
- 6DF - 3F 2 ) -j{D 2 + 3F 2 ) q jt - (D - F) 2 q r 
-l(D 2 + 3F 2 ) (l + 3q r ) 
\(D-F) 2 {2 + 2,q 3l ) 
l(D-F) 2 (l + 3q r ) 
1(D 2 + 3F 2 ) (2 + 3^-,) 
| {D 2 + 3F 2 ) (1 + 3g r ) 


(~ 27 + 9^0 

Hi + ll^r) C 2 
(-27 ~ 5*0 ° 2 

(~27 _ 9*) ^ 2 
(~27 ~~ 18*') ^ 2 
( — 5i ~ TElr) C 2 


~l + 
-| + 2q jt + g r 

| + 2g r 
' 2 

-3 ~ ??' 

- 3 — q r 

— 3 — 2q r 



TABLE IX: The coefficients P x , P' x , and A x in SU{3) flavor PQxPT for the A. 



Px P'x A x 



2 f {7D 2 -12DF + VF 2 )(l-3q ]l ) (_^ + i^) C 2 -§ + 2<&, 

^ (5D 2 + 30DF - 9F 2 ) -\{D + 3F) 2 gjl - § (7D 2 - Y2DF + 9F 2 ) q r + ±g r ) C 2 £ + g j7 + 2g r 

-±(D + 3F) 2 (l + 3cir) i + o. 



^i/> • 3F) 2 (1 • 3q r , ., , .,, 

- 2T ( 7Z)2 - 12Z)i? + 9i?2 ) (1 - %<) (ii - c2 I " 2 ^ 

-^(7D 2 -12 J DF + 9F 2 )(l-6g r ) (|4« r )C 2 \ - 2q r 

±(D + 3F) 2 (2 + 3 qjl ) -I -Qjl 

±{D + 3F) 2 (l + 3q r ) -3-* 



QxPT to do this extrapolation. QQCD, however, has no known connection to QCD and is 
of merely academic interest. Observables calculated in QQCD are often found to be more 
divergent in the chiral limit than those in QCD. This behavior is due to new operators 
included in the QQCD Lagrangian, which are non-existent in QCD. For the octet meson 
and baryon charge radii we find that such operators enter at NNLO. Hence, formally our 
NLO result is not more divergent than its QCD counterpart. This, however, does not mean 
that our result is free of quenching artifacts. While the expansions about the chiral limit 
for QCD and QQCD are formally similar, < r 2 >~ a + f3\ogrriQ + . . . , the QQCD result 
is anything but free of quenched oddities: for certain baryons, S~ and E~ in particular, 
diagrams that have bosonic or fermionic mesons running in loops completely cancel so that 
(3 = 0. In other words, < r 2 >~ a + . . . and the result is actually independent of mg! 
The same behavior is found for the charge radii of all mesons in QQCD as the meson loop 



TABLE X: The coefficients X , P' x , and A x in SU(3) flavor PQxPT for the AS transition. 



X 


Px 


P'x 


A x 


TT 

K 

uj 
ur 


-in D2 

^D(D - 3F) 

ir^ D ~ 3i? ) 


1 n2 

1 n2 
12a/3 
1 n2 

1 n2 
12-V/3 
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TABLE XI: The coefficients (3% and in 577(3) flavor QxPT for the octet baryons. 





8 Q 
Px 

IT 


K 


Px 

IT 


K 


n 

T) 

/ b 


3^Q 


n 
n 

u 


l/>2 
6 L Q 


n 

n 
u 


E+ 





_4 n2 





l/>2 


E° 
E" 






2 n2 
Z U Q 








1 />2 
12 U Q 




A 






2 n2 
3^Q 








1 n2 
12 U Q 




S° 





4 n2 
3-^Q 





1/72 

6^ 


EA 




^_n2 

3v^Q 


1 n2 


1 /72 

12^°<9 



contributions entirely cancel. 

PQQCD, on the other hand, is free of such eccentric behavior. The formal behavior of 
the charge radius in the chiral limit has the same form as in QCD. Moreover, there is a 
well-defined connection to QCD and one can reliably extrapolate lattice results down to 
the quark masses of reality. The low-energy constants appearing in PQQCD are the same 
as those in QCD and by fitting them in PQ^PT one can make predictions for QCD. Our 
PQxPT result will enable the proper extrapolation of PQQCD lattice simulations of the 
charge radii and we hope it encourages such simulations in the future. 
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APPENDIX: CHARGE RADII IN SU{2) FLAVOR WITH NON-DEGENERATE 
QUARKS 

In this section, we consider the case of SU(2) flavor and calculate charge radii for the 
pions and nucleons. We keep the up and down quark masses non-degenerate and similarly for 
the sea-quarks. Thus the quark mass matrix reads rn S Q^ = diag(m u ,md,mj,mi,m u ,md). 
Defining ghost and sea quark charges is constrained only by the restriction that QCD be 
recovered in the limit of appropriately degenerate quark masses. Thus the most general 
form of the charge matrix is 

Q su ^ = diag Q, ~, Qj , q h Qj , q^j . (A.l) 

The symmetry breaking pattern is assumed to be SU(A\2) L <g> SU(A\2) R <g> U(l)v — ► 
SU(A\2) V <g> U(l) v . 
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For the 7r + , 7r , and 7r° we find 



G P + V) = 1+ ' 



16tt 2 / 2 



g + ft J Fdd + ( g - gj J i^w - (1 - Qj + «) F ud - (- + qj j /•;,/ 



l + qi )F ld -( 2 -- q AF JU -( 2 -- qi 



lu 



+ a 9 —q 2 , (A.2) 



g pq = -G P J, and = 0, respectively. 

The baryon field assignments are analogous to the case of SU(3) flavor. The nucleons 
are embedded as 

B ljk = -j= (eijN k + e lk Nj) , (A.3) 

where the indices i,j and k are restricted to 1 or 2 and the SU{2) nucleon doublet is defined 

as 

"=(') {AA) 

The decuplet field %j k , which is totally symmetric, is normalized to contain the A resonances 
Tijk = %jk with i, j, k restricted to 1 or 2. The spin-3/2 baryon quartet is then contained 

as 

T m = A++, T 112 = i=A + , T 122 = i=A , andT 222 = A-. (A.5) 

The construction of the octet and decuplet baryons containing one sea or one ghost quark 
is analogous to the SU(3) flavor case |3j| and we will not repeat it here. 

The free Lagrangian for B and T is the one in Eq. (J2(Jj) (with the parameters having 
different numerical values than in the SU(3) case). The connection to QCD is detailed 
in (3j| • Similarly, the Lagrangian describing the interaction of the B and T with the pseudo- 
Goldstone bosons is the one in Eq. (j20j) . Matching it to the familiar one in QCD (by 
restricting the Bij k and %j k to the qqq sector) 

C = 2g A NS»A ll N + g 1 NS^Ntr(A fjl ) + g AN {ff A^Nje kl + h.c.) (A.6) 

one finds at tree-level 

4 1 „ 2 1 , „ . 

« = ^9a + -9i, P = -9i ~ -9a, and C = 9 an- (A.7) 

The contribution at leading order to the charge radii from the Pauli form factor F 2 (g 2 ), 
involves only the magnetic moments which arise from the PQQCD Lagrangian |3J| 

C = [ai« (B[S„ S V ]BQ SU ^) + n (B[S„ S U ]Q SU ^B) 

+ /I, str(Q su ^) (BIS,, S U ]B)} F>». (A.8) 

Note that in the case of SU(2) flavor the charge matrix Q is not supertraceless and hence 
there appears a third operator. In QCD, the corresponding Lagrange density is convention- 
ally written in terms of isoscalar and isovector couplings 

if 

C= 2M^ 0"°^' S ^ N + ^ N l S "> ^]r 3 iV) (A.9) 
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TABLE XII: The coefficients f3 x , 0' x , and A x in SU(2) flavor PQxPT for the proton. 



X 


0x 


Px 


A x 


uu 
ud 
dd 
uj 
ul 

dj 
dl 


l(Ag 2 A + 2g 1 g A +g 2 )(2-3q J ) 
-|<7a(5 + 6gz) + l>gi9A{2 - 3 Ql ) + \gl (1 - 9^ - 6 W ) 
-§ 5 ?(l + 3<ft) 
-l(Ag 2 A + 2g 1 g A +gj)(2-3q j ) 
-l(4g 2 A + 2g 1 g A +g 2 1 )(2-3q l ) 

yui+sq 3 ) 

ig 2 (l + Sqi) 


(~ 27 + TElj)9AN 
(-§ + 5<7j + TSlfiffAN 
(27 + 9<1i)9an 
(27 ~~ j8lj)dAN 
(27 ~~ 1s1i)9an 
( _ 27 — 9Qj)9AN 
(~27 ~~ 9®) 9 AN 


-l + 2 qj 
Qj + 2qi 

l + m 
1 - 2g, 

-5-?? 
~3-« 



and one finds that the QCD and PQQCD coefficients are related by {31] 

A*o = 7: (A^ + + 2// 7 ) , and //1 = - (2/io, - //^g) . (A. 10) 

o 

Likewise, the leading tree-level corrections to the charge-radii come from the Lagrangian 

£ = [c Q (Bi3Q sc/ ( 2 )) + c p (BQ SU{2) B) + c 7 str (Q SU ^)(BB)] v^F^ (A.ll) 

x 

that matches onto the QCD Lagrangian 



C = T2 ( c ° NN + c i Nt3N ) v A F ^ ( A - 12 ) 

X 

with ^ 1 

Co = - (c„ + + 2c 7 ) , and ci = - (2c a - cp) ■ (A. 13) 



Evaluating the charge radii at NLO order in the chiral expansion yields 



2 6c 3a 1 ^— \ 

^ T E aT + o T\ ,r2 1 a„2 f 2 / -> 



A2 2M& 16tt 2 / 2 



2 2 

A x log ^ - 5 log + 10 fte(m X) A, //) 



(A.14) 

The coefficients c are given by c p = Co + ci and c n = c — ci while a p = /io + /ii and 
a„ = fiQ — fii. The remaining coefficients are listed in Table IXTTl for the proton and Table IXTTTl 
for the neutron. 
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